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Abstract: We compute the anomalous dimension for a short single-trace operator in 
planar ABJM theory at intermediate coupling. This is done by solving numerically the 
set of Thermodynamic Bethe Ansatz equations which are expected to describe the exact 
spectrum of the theory. We implement a truncation method which significantly reduces 
the number of integral equations to be solved and improves numerical efficiency. Results 
are obtained for a range of 't Hooft coupling A corresponding to < h{\) < 1, where h{X) 
is the interpolating function of the AdS4/CFT3 Bethe equations. 
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1. Introduction 

Recently there has been a significant interest towards integrable structures which arise 
in the context of AdS/CFT correspondence with the best-studied example being the 
AdS5/CFT4 duality between four-dimensional planar = 4 Super-Yang-Mills theory 
and Type IIB superstring theory on AdSs x S^. There are also other AdS/CFT dual 
pairs [Q, ^ where integrability gives us serious hopes of solving exactly a highly non- 
trivial quantum field theory. One of these is the ABJM duality proposed in which 
relates three-dimensional planar J\f = 6 super Chern-Simons theory and Type IIA string 
theory on AdS4 X CP^. It appears that the spectrum of anomalous dimensions / string state 
energies in this AdS4/CFT3 duality may be found exactly at any value of the coupling 
m by applying an approach similar to the one used in AdS5/CFT4. In particular, the 
Bethe ansatz equations, which describe the spectrum for asymptotically long single-trace 
operators at any coupling, have been proposed in |p and extended to all loops in Q, 
. These equations however do not capture the so-called wrapping interactions Q which 
means that other tools have to be used in order to obtain the spectrum for short operators 
or the energies in finite volume. 

In the AdS5/CFT4 case this issue has been successfully addressed by means of the 
generalized Luscher formulae |^ and the Y-system/Thermodynamic Bethe ansatz (TBA) 
approaches |jl^, [0]. For AdS4/CFT3 a Y-system of functional equations was proposed in 
||To|] , and was later refined as well as supplemented with a set of TBA integral equations in 
1 12], [^]. Unlike the Bethe ansatz, the TBA and Y-system are expected to be valid for any 



state, thus providing a way to obtain, in principle, the full exact spectrum of the theory. 
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For AdS5/CFT4 the Y-system and TBA have passed a number of nontrivial tests. 
In particular, the Y-system allows one to efficiently reproduce perturbative gauge theory 
calculations at weak coupling (see e.g. fl^), while at strong coupling it matches the 
semiclassical predictions obtained from the algebraic curve [14|, [15|. In addition to these 
analytical checks, important numerical results have been obtained from the TBA in |16|, 
1 17] where the anomalous dimension of the Konishi operator was computed for a wide range 
of values of the 't Hooft coupling. The strong-coupling predictions obtained in these works 



seem to agree with the analytical results obtained by several other methods |18], ||l5| (see 
also |19|), providing yet another successful test of the proposed TBA and Y-system. Also, 
very recently the TBA equations were reduced to a finite set of integral equations [20|. 
Analogous checks have been done for AdS4/CFT3 - the four-loop wrapping corrections 



obtained in |1C] were reproduced in |^], [Q, while in [13| the proposed Y-system and 
TBA were shown to be remarkably consistent with the algebraic curve quantization |23|. 
However, a numerical analysis similar to ||l^, |17] has not been attempted so far and would 
be an important test of the integrability properties in this theory. 

Here we present a first step in this direction, solving numerically the TBA equations 
for one of the simplest unprotected operators in the s/(2) sector to compute its anomalous 
dimension non-perturbatively. We start from weak coupling and are able to reach those 
values of the 't Hooft coupling A for which the AdS4/CFT3 interpolating function h{X) 
becomes equal to 1 (this should correspond to A ~ 1 as well) ^. To facilitate efficient 
numerical analysis, we implement a truncation method [24| for the TBA equations which 
is based on partially solving the underlying Y-system and allows us to eliminate from the 
equations an infinite number of unknown functions. 

This paper is organized as follows. In section 2 we introduce the state that we are 
studying, in section 3 we present the TBA equations, and in section 4 describe the trun- 
cation method. In section 5 we discuss the numerical results, and we conclude in section 
6. Appendix A contains notation that we use. 



2. Description of the state 

The gauge theory operator that we study in this paper is the L = 2 state in irrep 20 
of SU{4:), described in the Appendix of Grouping the scalar fields of the theory into 
SU{4:) multiplets as follows: 

y^ = {Ai,A2,bI,bI) y] = (4,4,i?i,i?2), (2.1) 

we can write this operator as 

o = tr (y^yiy^y])x^g, (2.2) 

where the coefficients Xcd antisymmetric in both A, B and C, D pairs of indices (for 
more details see This is one of the simplest unprotected operators in the Af = 6 

^Since the TBA equations include h{\) rather than A, our result is the anomalous dimension as a function 
of h{X) and not of A. 
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super symmetric Chern-Simons theory. In the asymptotic Bethe ansatz (ABA) of this 
state is described in su{2) grading by two momentum-carrying Bethe roots - one M4 and 
one U4, which are equal. The Bethe equations in this case reduce to^ 



^ =1 (2.3) 



where the function x{u) is defined as 



with standard two branches 




-*(") = H^VS-'V^') . -"» = il? + n/4-^.|. (2.5) 

and in the Bethe ansatz equations we use the physical branch. We also used the general 
notation 

/± ^ fiu ± i/2), ^ /(n + m/2) . (2.6) 

The function h{X) in ( |2.4| ) is the so-called interpolating function (see |2^) which plays the 
role of the effective coupling in the Bethe ansatz and TBA. Its weak and strong coupling 
expansions are known to be 

h{X) = A + hsX^ + O (A^) = 7X72 + h^ + (^-^ . (2.7) 

The weak-coupling coefficient /13 = —8 + 2({2) was computed directly from perturbation 



theory in [^], ||2^. For the strong-coupling coefficient h several calculations suggest 



different values: ||2^, |Q, [^] argue that it is zero, while propose the value 

-i^ (see also ||)3. 

The TBA equations correspond to the sl{2) rather than su{2) Bethe ansatz equations, 
so we need to make a duality transformation in the Bethe ansatz. We find that the 
sl{2) representative of this state is described by the same two Bethe roots but has L = 1 
rather than L = 2, the corresponding Bethe equation being 

L 

= -1 . (2.8) 



1-4 

The corresponding single-trace operator is in the same supermultiplet as ( |2.2D (and has the 
same anomalous dimension). Its explicit form can be found in |3^] ^ and its bare dimension 
is 3. 

Importantly, from the ABA equations (2.8) we see that for any coupling the two roots 



remain exactly at zero^: ii4 = = 0. This means that, unlike what happens for Konishi 



^the dressing factor a does not appear because a{u,u) = 1. 

^Related calculations at strong coupling on the string side of the duality have been done in [ ^ . 
^Roughly speaking, it is a linear combination of terms with two scalar fields and two covariant derivatives 

^This is also consistent with the zero-momentum constraint Q: Y[ 11 -^ = 1- 



-3- 



in the = 4 SYM case, the Bethe ansatz result for the scahng dimension can be found 
exactly (in terms of the interpolating function h{X)). It is given by 



Eaba = Eo + ef{u4) + ef (nj) (2.9) 
where Eq is the bare dimension of the operator and 

- hiX) (-^^ - -j^) , (2.10) 

which means that in our case 



Eaba = \/lQh{Xy + 1 + 2 . (2.11) 

It is well-known that the ABA result is usually incomplete because of the wrapping 
interactions which arise due to the finite length of the corresponding spin-chain Q. The 
exact anomalous dimension includes a correction, 6E, to the ABA result: 

E = Eaba + SE . (2.12) 

The leading weak-coupling term of this correction (first wrapping) has been computed in 
p!0| ] and confirmed in ||2^ (see also [36] where wrapping corrections for similar operators 



were studied). It has the form^ 



SE^rapping = (32 - 16C(2))/i(A)^ + O(A^) (2.13) 



-'wrapping 

and thus 

Eaba + SE^rapping = 3 + 8h{Xf - 16C(2)/i(A)4 + OiX^). (2.14) 

Our main result in this work is a numerical computation of the correction 5E for < 
h{X) < 1 by means of the Thermodynamic Bethe Ansatz. 

3. TBA equations for AdS4/CFT3 

The Y-system which describes the spectrum of the ABJM theory in the planar limit was 
first proposed in ||l^ and later refined in [ 13 1 , ]l2| . This system of functional equations can 



be summarized in the diagram shown in figure |^. 



We denote the various Y- functions as in |13| - there are fermionic functions 1© and Ig 



bosonic functions Yq^ (n = 2, 3, . . . ), pyramid functions Ya„ (n = 2, 3, . . . ) and middle node 
functions (n = 1, 2, . . . ). As the state we consider belongs to the sl{2) subsector of 

the theory, the middle node Y-functions of two types, corresponding to two series of black 
nodes in figure |l|, are pairwise equal. We denote these Y-functions by or Ya^, so that 
we have = Y^ = = Yafi for all a. 



®since at weak coupling h{\) = A + 0(A^) this expression can be rewritten also as SE = (32- 16C(2))A'* + 
0(A5) 
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Figure 1: A graphical representation of the Y-system proposed in [Q for ABJM theory. Each 
circle in this infinite 3D lattice corresponds to a Y-function. 

The TBA equations which describe the exact energy of the ground state in finite volume 



were proposed in [|^, [13|. In [|13[ they were also extended via the contour deformation 
trick [p^ to excited states in the sl{2) subsector, the resulting equations being: 



logF« = +Km^i * log + 27^^°^) * log(l + li,J + 2 



logl® = -Krn-1 * log 



1 + 

1 + yYo„ 



log 



R 



{+) 



2<^)*log(l + F^J-2 



1 + 

loglA„ = -Kn-l,m-l * log(l + Ya^) - Kn-1 * lo, 



log- 



B 



(+) 



B 



(-) 



1 + 1/1® 



+ i7r (3.1) 

- ivr (3.2) 
(3.3) 



+ 2(7^(°^)+^fM*log(l + y^) 



+ 2 



2 



v-^ , R[^\u + ik) 


n-i Rl {u + ik) 



+ 2 



n — 3 
2 



^ , Bi-^\u + ik) 

2^ lor 



logyo„ = Kn — l,m— 1 *log(l + l/yo„)+i^n-l*log 



i + y» 



logy^ = (L + K4)iog- 



1 + i/y® 
* iog(i + 1/Yb) + 7^i\°^ * iog(i + 



(3.4) 
(3.5) 



+ (7^(l^)+<2-2)*log(l + yA^) 

+ f25„^ - n^^i^ + aii^)) * log(l + Y.J + 



71-1 

2 



log <I>(u + i/c) 



2~ 



+ ivrn 



Here and throughout the paper we use notation that is given in Appendix A. We also 
assume summation over the repeated index m. 

The exact Bethe roots u^j = u^ j are fixed by the exact Bethe equations, 

y^^n4,,) = -1, j = l,2,...,i^4 (3.6) 
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and in general the values of these Bethe roots may differ from those one gets from the 
ABA. In the case we are studying the roots remain at zero within our precision - we have 
checked this numerically, verifying (3.6) with the use of equation (^]^) which we analytically 
continued as in |16|. 

The energy of the state we are studying is written in terms of the Y-functions as 

pGO ^ o mir z' \ 

E = E.«. + .E, .E = 2X:/^|ii?!^Mi + y-) ,3.7) 



where ea{u) is defined by ( ^.10 ) 



As our goal is to solve the TBA equations numerically, we will make to them several 
modifications, which are described below. 

First, we substract from the original equations the equations which are satisfied by 
the asymptotic large L solution of the Y-system [0], Let us describe this solution 

for the state discussed in section 2 which belongs to the sl{2) subsector and has L = 1, 
K4, = Ki = 1 with two Bethe roots ti4 = nj = 0. We use bold font to denote the asymptotic 
Y- and T-functions. The main formulas are: 

Ya. = ^ "-^"'^ - 1 , 1/Yo. = ^ - 1 (3.8) 

^ _ T2,3T2,1 _ T\2Ti£ 

-'-3,l-'-l,2 -'-2,0-'-0,l 
All the 'Ta,s functions can be found from the generating functional W. In particular 

00 00 
W = ^Ti,,(u + ii^)D^ , W-^ =Y,{-'^TT^aAu + i^)D^ . (3.10) 

s=0 a=0 




-1 



where D = e . Finally, 



(3.11) 



^•M= I ^ I Mu)Ta,i{^). (3.12) 




A remarkable property of these asymptotic Y-functions is that they satisfy a set of 
TBA-type equations. This is true not only for A — )• 0, but also for finite values of A and L 
(see e.g. [pO], p8|)'^. The equations satisfied by the asymptotic solution are the same as 



''For the state we are studying we have also checked this numerically 
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the original TBA equations given above, except that the terms in the r.h.s. which involve 
the li^ -functions should be omitted. After substracting these equations from the original 
ones and also slightly rewriting the kernels for more convenient numerics (similarly to ) 
we get: 

log ^ = +Km-i * log V'r°"/J"Jv"' + * log(l + Y^) (3.13) 

log ^ = * log - 2B^^"^^ * log(l + Y^) (3.14) 

I© i + Ya^ -1 + 1/ I o„ 

log = -Kn-l,m-l * log - Kn-1 * log , (3.15) 

Ya„ (1 + Ya,„) 1 + 1/Yb 1 + 



+ 2Mnm * log(l + Yn 



^ = * (1 + 1/Yc.) + ^"-^ * 1 + 1/ye 1 + Y. ^'-''^ 

log^ 



+ J\fnm * log 



(1+Ya„) 

+ (2Snm - n^^^ + B'i^^) * log(l + Y^] 



Note that in the Konishi case the function 1a2('u) had poles at the positions of the Bethe 
roots; in accordance with that, in our case it has a double pole at n = 0. The asymptotic 
Y-function Y/s,^{u) also has a double pole at n = 0, so the combination l^^^^ ^^^-^ which 



appears in the equations ( 3.13 )-( 3.17 ) has no singularities on the real axis 



In the next section we discuss a further simplification to these equations which leaves 
only a finite number of unknown functions. 

4. Truncating the TBA equations 

The truncation method which we describe in this section has been proposed in and is a 



simpler version of the general treatment in |2^. Unlike the method of [20| it involves some 
approximations, and it also relies on certain analyticity assumptions for the Y-functions. 
We will not discuss these assumptions here, but we have confirmed numerically that in our 
case the resulting equations are consistent with the original TBA of []T3| . 

The truncation is done as follows. First, remarkably, it is possible to eliminate all the 
Yq^ functions from the TBA equations, replacing them by a single unknown function. The 
reason is that the Y-system equations for these functions are quite simple, 

^ci^d: = (l + Yo„+J(l + Yo„.J , n>3 (4.1) 
and to solve them we can use the following ansatz for the corresponding T-functions (see 

Ti,s = s + Ks*fR (4.2) 

*see also H 
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where /r is a new unknown function. The lo^ functions are obtained |1C], |13], as usual, 
from ^ 

Yo. = ( /''''l'''' -l] , (4.3) 



and it can be shown that the infinite set of equations (4.1) is indeed satisfied for any //j 



with suitable analyticity properties. This function is fixed by the TBA equations - the 
infinite set of equations ( |3.16D reduces to a single equation for That equation has the 
form 

1 + y» _ (1 + A'+ fn + fR/2) (1 + Jr + fn/2) 



1 + 1/1® (1 + K+ fn - fR/2) (1 + K- Jr - fR/2) 



(4.4) 



where *p.^. denotes principal value integration. Solving for fR we get an expression suitable 
for numerical iterations^: 

^ + l)(;/ + 1) - V4(a + 1)2^2 _ 52(^ _ 1)2 ^^^^^ 

2/ - 1 

where y = j^^, o = Re (/C^ *p.^. /a) , 5 = Im (/C^ *p.„. /r). 
Note that the Y-system equation for Yq^, which reads 

y^yo,- ^^yy^ , (4.5) 

is not satisfied automatically for arbitrary fR (since it includes the fermionic functions 
1®, 1®), but will hold provided that the TBA equations are satisfied. 

So far we haven't made any approximations in the TBA equations - the truncation we 
have just described is exact. However, we cannot directly apply the same method in the 
upper wing of the Y-system to replace the functions X\„,Ym^ by a finite number of unknown 



ones, as the Y-system equations for 1a,j,1«„ are more complicated than (4^). But since 
the li^ functions decay fast with n, within our precision it is enough to keep only the first 
6 or 7 of them and set the other ones to zero. With this approximation we see that for 
n > M, where M the number of li^ functions that we retain, the pyramid Y-functions Ya„ 
decouple from the rest of the Y-system and are governed by an equation very similar to 



This equation is solved by an ansatz analogous to (4.2): 



Ta,l = a + Ka-M+l * fu (4.8) 



with Ya = 7r — "''^^"'"^ 1. Here fu is another new unknown function, and it is fixed by 

an equation similar to (|4.5| ): 



fu = a(4y - 2) - 2{^ {y - l)y{2a + M + 1)^ - - {M + l)y + M) (4.9) 



'Note that fu is only nonzero for —2h <u< 2h. 
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where y = 1 + Ia^ , a = Re {K^^ *p,y, fu), b = Im {K^ *p,^, fu) ■ 

Lastly, we need to rewrite the r.h.s. of the remaining TBA equations in terms of the 
new functions fR,fu- To do this, we wih need the exact T- functions - which are given by 



|4.2| ), (4.8) - and the asymptotic T-functions, which are obtained from the same expressions 



when fn, fu are replaced by f/j, 

W = ^mir+ _ ^^^mir- ' ^u{u) = Km-i{u,v)* (l/^mir+^ _ ^mir- • 

(4.10) 

We list the remaining TBA equations below. 
Equations for fermions: 

log ^ = log Ki* log Km-1 * log — + Km * log — (4.11) 

I® -'-1,2 i-M+1,1 -Lm,1 

1 + 



K^m-i * log Y^Y^ + 2^^°"^ * + ^•^^ 



m=2 



log — = - log + Ki * log — ^ + Km-1 * log = ■ Km * log 

I® -11,2 -Ll,! -l-M+l,! i-M,l 



1 + Y. 



+ ^ i^„_i * log - 2^(0-) * log(l + Y^). (4.12) 

m=2 ^"^ 

Equations for pyramids (n = 2, . . . , M) 

1 s- 1 Ta/+i 1 T/vz+i 1 Tjvf,! 

log = -0„,A/ log — Kn-l,M-l * log — h Kn-1,M * log =^ 

IA„ -l-AZ+l,! i-M+1,1 i-M,l 

- 2^ ^n-l,m-l * log — — p Kn-1 * log . , . ,y -, , ^ (4-13) 

+ 2A1„„*log(l + y,„). 
Equations for middle nodes (n = 1, 2, . . . , M): 

log ^ = f25„^ - 7^1l^) + ei^^)) * log(l + y^) (4.14) 



+ 7^("0)® f Km-1 * log 1^ - Km * log ^ + £ * log ^^^] 

+ K-l,M-l * - ^n-l,M * log TfT^ + 2^ ^r-l,m-l * log , ^ 



As a result, the equations that we are solving numerically are ([4.5|), ( ^.91) and ( 4.11] )-( ^.14| ). 
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Figure 2: The correction SE to the Bethe ansatz resuh for the energy, as a function of h{X), is 
shown by black dots. The solid line is the first wrapping correction given by ( 2.13| ). 



5. Numerical results 



Our numerical strategy is similar to [| 
Q, (pID -( ^) by iterations'^ 



and we are solving the TBA equations ([4.5|), 
We truncated the number of middle node Y-functions 



to M = 6 or 7, and used cutoffs on the u axis to make finite the range of integration in 
the convolutions. The range of coupling that we study is < /i(A) < 1 (since the TBA 
equations are written in terms of h{X), we obtain the anomalous dimension as a function 
of h{X)). 

Our main result is the correction 6E to the ABA 
value for the energy (see ( 2.12| )) and it is shown in 
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Fig. |2[ The numerical values of this correction are 
also listed in table |^. We estimated their absolute 
precision as about ±10"^ for h{X) < 0.55 and ±10"^ 
for h{X) > 0.55. In Figure |3| we plot the full energy, 
including the ABA part. As expected, at weak cou- 
pling our numerical results are completely consistent 
with the leading wrapping correction obtained ana- 
lytically in [|o|, |2^ (this is also seen in Figure |2|). 
As the coupling is increased, our results gradually 
deviate more and more from that prediction. 

The behaviour of middle-node Y-functions 1«„ 
exhibits several interesting features. While in the 
AdS5/CFT4 Konishi case all these functions take pos- 
itive values, here they have signs alternating with n, 
i.e. l«j < 0, li^ > 0, 1«3 < etc. In Figure § (left) we plot their absolute value for 
h{X) = 0.9. Since these Y-functions appear as log(l + Ym^) in the expression for the energy 
and in the TBA equations, their negative values beyond —1 would give rise to singularities. 

^"Sometimes we also had to iterate the equation for 14^ separately from the others to make the process 
converge. 



Table 1: Numerical values of the cor- 
rection to the Bethe ansatz result for 
the energy, with uncertainty in the last 
digit. 
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Figure 3: The full scaling dimension Eaba + 5E. Dots: numerics, solid line: ABA, dashed line: 
ABA + 1st wrapping as given by ( ^.14 ) 
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Figure 4: Left: The middle-node Y-functions Yn.o{u) for h{\) = 0.9. The figure shows plots of 
minus Yi^ (blue), ¥2,0 (purple), minus ^3,0 (yellow), 14^0 (green) and minus 15,0 (red). Right: the 
exact Y-function Yix> (solid line) and its asymptotic counterpart Yi^o (dashed line) for h{X) — 0.9. 



We discovered that as the coupling is increased, the middle-node Y-function indeed ap- 
proaches the critical value —1 for u close to zero (while other Y-functions remain far from 
this dangerous value). That is clearly seen from Figure ^ (left), and at h{\) = 1 we have 
Y%^{Q) ~ —0.993 which is already very close to the singularity. This means that proceeding 
further in the coupling would probably require obtaining with a high precision. 

The corresponding asymptotic Y-function already starts to take values beyond — 1 
at /i(A) ~ 0.6. However, that does not lead to any singularity since log(l + Y«j) never 
appears in the TEA equations. In Figure ^ (right) we plot for comparison the exact and 
the asymptotic Y-function. 

It is possible that if the coupling is increased further the function will cross the 
critical value —1. Then the TEA equations will probably require some modification such 
as extra driving terms, and it would be very interesting to understand whether this indeed 
happens for the state we are studying. For other models similar issues have been explored 
in 1 37], while in the AdS/CFT case the possibility of such singularities arising was discussed 



in e.g. ig, g, [||. 

The second middle-node Y-function Y^ also shows unusual behaviour, rapidly increas- 
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Figure 5: The middle-node Y-functions Yi_o(u) (left) and l2,o(u) (right) for several values of the 
coupling: h{X) = 0.6 (blue), 0.7 (purple), 0.8 (yellow) and 0.9 (green). 





Figure 6: The functions fR{u) (left) and fu{u) (right) for h{X) — 0.9 and M — 7, shown by solid 



lines. The dashed line in the left figure shows the asymptotic ffl function (4.10), while fu almost 
coincides with its asymptotic expression fjy. 

ing in magnitude as the coupling is being increased. This is shown in Figure ^ (right). 

Lastly, in Figure ^ we show the plots of the new functions //j, fu which parametrize 
the solution of the Y-system in the right wing and in the upper wing, respectively. 

The total range of the coupling we have investigated, < h{X) < 1, is four times greater 
than the convergence radius |/i(A)| < 0.25 of the weak-coupling expansion of the ABA result 



( 2.11 ), which suggests that we are exploring the intermediate coupling regime. We were not 
able to make a consistent prediction for the strong coupling expansion coefficients, but we 
hope this could be done in the future by going to greater values of the coupling. Increasing 
the coupling poses a challenge because beyond the value /i(A) = 1 the iterative procedure 
we use for solving the TBA equations converges too slowly; hopefully this problem may be 
overcome by improving the numerical algorithm (e.g. using Newton's method). 

At strong coupling the leading term in the exact energy should be proportional to 
A^/^ II which in our case is equivalent to \/h{\) (while string theory predictions for 
subleading terms are not available as of now). However, the strong-coupling behavior of 



the ABA result ( 2.11 ) is completely different 



Eaba = 2V2X + 0{1) =Ah{X) + 0{l), (5.1) 

suggesting that the leading asymptotics predicted by the ABA should cancel against^^ the 
^^In contrast, in the AdS5/CFT4 Konishi case ]l6| , [l7| , |l8| , [iHj the ABA result already has the correct 



- 12 - 




Figure 7: A log-log plot of the full conformal dimension E vs. h{X). We show the ABA prediction 
(soHd curve), its asymptote (dotted line defined hy E = 4/i(A)), our numerical data (black dots) 
and the expected slope of the result at strong coupling (dashed line defined hy E = const x y^h{X)) . 



TBA correction 5E. It is possible that we can already see this starting to happen for 
/i(A) > 0.8 when the exact energy becomes smaller than the ABA prediction. In Figure 
1^ we show a log-log plot of the full energy, and the ABA result asymptotes to a straight 
line, consistently with (^.1|), while one may expect that the exact energy will asymptote to 
a straight line with a different slope^^. 



6. Conclusions 



In this paper we have applied the Thermodynamic Bethe ansatz approach to compute the 
exact anomalous dimension of a short operator in AdS4/CFT3, for the first time solving 
numerically the TBA equations for this theory. We have explored the values of the coupling 
< h{X) < 1 and at weak coupling our results are consistent with known predictions. We 
also found that as the coupling is being increased one of the Y-functions approaches the 
critical value —1. It would be very interesting to go to higher values of the coupling, 
and achieving that may be possible by improving the numerical strategy or perhaps by 
applying the very recent approach of [20|. Hopefully a comparison with string theory 
calculations could be made eventually at strong coupling, providing a deep test of the 
integrable structures in AdS4/CFT3 correspondence. 
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7. Appendix A: Notation 



The R and B functions are 



= n H«) - 



^i*' = n 



x{u) 



(7.1) 



where for the state we consider in this paper all Ka are zero except K4 = = 1. The 
scalar factor $ is defined as 



$(u) 



5(+)+^(-)- ( ^^+\ ^4 



1-1 

2 



n^'(^'^4j), M^)= n H^+ik). (7.2) 



The kernels in TEA equations are: 



. . 1 d ^ u - V + in/2 
Kn{u,v) = — ^In 



2ra/7r 



2'!ri dv u — V — in/2 n^ + 4:{u — v)"^ ' 

m— 1 n— 1 
2 2 

^n,m(^i,^^) = XI XI -f^2j+2fe+2(^i,^^), 



_ m—1 I, n— 1 

m — 1 n — l 

2 2 



(7.3) 
(7.4) 

(7.5) 



m—1 I, n—l 
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where we assume -fCo,n = 0, K^^ = 0. Also, 

(7.6) 

Snm{u,v) = -L|.ioga5^5(xI+"](n),x[-"lH,x[+'"](^;),x[-'"](t;)), (7.7) 

Sn„^{u,v) = Snm{u,v) + '^V^'^\v), (7.8) 



where 



and 



n — l m— 1 

2 2 



B^-^)(uv)= V V J- A w + + ^ ~ ^^/^ + (79) 

3— 2~ ^~ 2~ 



n — l m— 1 

2 2 



^-^ ^-^ ^ l-Ki ov r{u — la/ 2 + ij,v + 10/ 2 + ik) 



2 2 



o [+n] _ [— m] -, 

7^M(n,t,) = — log ^r 1 1 --p('")(^^), (7.11) 

A^nm = i^„_l(fe7^(°") + , (7.13) 

AA„^ = 7^("0)(|=)K„_l + , (7.14) 



-J-a„ log + (7.16) 

^ ' 2tt ^ x'^'^iv -ia/2) ^ ' 



The convolutions in integral equations are over the second variable, and * denotes 
standard convolution along the real axis: K{u,v) * f{v) = dvK{u,v)f{v). 

The symbol ® denotes convolution along with analytical continuation across the cut 
u e (-00, -2h) U (2/i, +oo) (see |3[, [||]). E.g. 

7^("°)®log(l + y«)= / dv 7^("0) iog(i + y«) - iog(i + i/y®) . (7.17) 

One should also be careful about which branch of x{u) to use in various places in TBA 
equations. The mirror branch is used when u is the free variable or the variable that is 
being integrated over. Otherwise the physical branch is used. 
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